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Abstract 

Spectrum of masses of light pseudoscalar (ir, K, r), rf) and vector 
(p, K*, u>, 4>) mesons can be explained on the base of the following 
assumptions: (1) analytical confinement (propagators of quarks and 
gluons are entire analytical functions of the Gaussian type), (2) mesons 
are bound states of quark and gluons (the Bethe-Salpeter equation in 
the one-gluon exchange approximation) and (3) QCD coupling con- 
stant a s (M) is a monotone decreasing function of mass of a bound 
state M. 

The decay constants and fx are calculated. 

1 Introduction 



At present time Quantum Chromodynamics (QCD) is considered as the true 
theory which control the behavior of quark and gluons as "bricks" of hadron 
matter. Due to the confinement phenomena in experiments we observe col- 
orless hadrons as bound states of quarks and gluons. Up to now a rigorous 
analytical solution of the problem of confinement and hadronization of quarks 
and gluons into hadrons is still missing. Therefore effective phenomenologi- 
cal and semi-phenomeno logical approaches based on QCD are developed to 
derive non-trivial statements about hadronic processes. Because dynami- 
cal mechanism of transformation of quarks and gluons into hadrons is not 
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clear now the conception of quark and gluon structures depends on phys- 
ical context under consideration and the aim of all theoretical approaches 
claiming to be obtained from "the first principles" of QCD is to find process- 
independent concepts using symmetry and group arguments (chiral, flavour, 
anomaly, mixing and so on) in order to get possibilities to compare differ- 
ent characteristics of physical processes (there is huge literature on these 
problems, see, for example, [1, 2]). 

We would like to stress that the application of methods of quantum field 
theory can be successful if (1) propagators of interacting particles are known 
and (2) the effective coupling constant is small. Namely this situation takes 
place in the standard model of electro-weak interactions. What have we in 
the QCD as theory of "strong interactions"? 

From physical point of view it is evident that processes of hadronization 
and confinement of quarks and gluons take place in the same space-time 
region, however the behavior of quarks and gluons in this confinement region 
is not known and, from our point of view, namely this circumstance is the 
main reason why QFT methods can not be applied directly. From our point of 
view if the behavior of quarks and gluons is known in the confinement region 
and the coupling constant is small enough, we will succeed in analytical 
solution of hadronization problem. 

The idea to get the quark propagator in the confinement region is not new. 
The Schwinger-Dyson equation is the main tool to calculate this propagator 
(see, for example, [3, 4, 5]). However the theoretical problem is what quark- 
gluon vertex and what gluon propagator should be used in this equation 
and then what computing methods should be applied to get the solution of 
this equation. As a result calculations become to be very cumbersome and 
opaque (try to repeat these calculations!). 

The confinement is generally accepted to be a result of nonperturbative 
nonlinear interactions of gluons in QCD. This idea plus the Wilson loop 
confinement argumentation are used to reduce the relativistic hadroniza- 
tion problem to a stationary Shrodinger picture with an increasing potential. 
Thus confinement is a static picture, i.e. there exists a constant in time po- 
tential which keep two quarks together. Conception of the confinement of 
one isolated quark is not formulated at all. 

Our point of view (see [7]) is that the cause of confinement is insta- 
bility of bosons in a homogeneous self-dual vacuum field. In particular in 
the quantum electrodynamics confinement can not take place because all 
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fundamental particles (leptons and baryons) are fermions. In QCD gluons 
(bosons!) owing to nonlinear self-interaction play double role - they trans- 
fer the interaction as photons in QED and play role of real particles-bosons 
with zero mass. Namely massless gluons being bosons leads to instability 
of QCD-vacuum which is realized for a nonzero homogeneous self-dual vac- 
uum gluon field which in turn leads to analytical confinement quark in this 
vacuum field. Thus from our point of view the initial "free" quark-gluon 
Lagrangian should contain this field and propagators of quarks and gluons 
described by this "free" Lagrangian should satisfy the confinement criterion. 
Physically analytical confinement means that quarks and gluons are fluctua- 
tions in space and time. Space-time scale of these fluctuations is defined by 
the strength of vacuum gluon field A. 

The second question, what is the QCD coupling constant in the confine- 
ment region, i.e. for low energies E < 1 Gev ? Usually coupling constant 
is supposed to be constant in this region although this statement is not 
completely consistent with general behavior of the running QCD coupling 
constant. Various investigations result in a remarkable variety of its infrared 
behavior (see [6]), so that this question can be considered to be open. Thus 
the unique answer this question is not yet exist, so that some speculations 
can be done. 

The main subject of this paper is the spectrum of masses of light pseu- 
doscalar (tt, K, rj, rj') and vector (p, K*, cu, (p) mesons. The aim of modern 
theoretical approaches is to get correlations between masses, i.e. so called 
mass formulas on the base of symmetry arguments (see, for example, [1]). 
Our point of view if the propagators of quarks and gluons are known in the 
confinement region and the QCD coupling constant a s (M) is small enough, 
then the Bethe-Salpeter equation can be used to calculate desired masses. 

The main result of this paper is that mass differences in pseudoscalar and 
vector multiples can be explained on the base of the following assumptions: 
(1) analytical confinement (propagators of quark and gluons are entire ana- 
lytical functions of the Gaussian type), (2) mesons are bound states of quark 
and gluons (the Bethe-Salpeter equation in the one-gluon exchange approx- 
imation) and (3) QCD coupling constant a s (M) is a monotone decreasing 
function of mass of a bound state M . 

Our approach is based on the following statements 

• The methods of Quantum Field Theory can be used, it means that 
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weak coupling constant regime should take place and perturbation cal- 
culations can be applied. 

• Our guess is that the selfdual homogeneous gluon fields with a con- 
stant strength is a good candidate to realize the vacuum QCD. In this 
fields analytical confinement takes place, i.e. propagators of quarks and 
gluons are entire analytical functions in the p 2 complex plane. 

• If propagators of constituent particles are known and coupling con- 
stants are small enough bound states can be found as solutions of the 
Bethe-Salpeter equation in the one-gluon exchange approximation. 

• The variation of the QCD coupling constant a s (M) in the low energy 
region M < 1 Gev should be taken into account. 

Our aim is to understand the general features of spectrum of light mesons 
in the most simple dynamical way so that we want to simplify the problem as 
much as possible. The first observation is that the quark and gluon propaga- 
tors can be approximated by virton fields, i.e. by pure Gaussian exponents. 
Besides solutions of the Bethe-Salpeter equation in this case can be found 
in the explicit analytical form, so that qualitative characteristics of the mass 
spectrum can be understood more profoundly (see [11]). The second ob- 
servation is that for light quarks located in a selfdual homogeneous gluon 
fields with a constant strength the main contribution into quark propagators 
comes from so called zero modes. We show that these two points defines 
main features of light meson spectrum. 

Thus our formulation of the problem looks: 

Does it exist a reasonable form of propagators of quarks and gluons in- 
duced by the behavior of constituents in a selfdual homogeneous gluon fields 
with a constant strength and the QCD coupling constant a s (M) in the region 
M n < M < M v i for which one can obtain the masses of pseudoscalar and 
vector mesons ? 
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2 Lagrangian and propagators. 

Our basic assumption is that the QCD vacuum is realized by a self-dual gluon 
field with constant strength 

B lt (x) = B°(x)1?, B;{x)=n a B ll (x), B^x) = A 2 b^x v , (1) 
Bp U (x) = d^B v {x) - d v B^(x) = -2 A 2 b^ = const. 

Here n a is a constant vector in color space. The parameter A defines the 
confinement scale and 

The field Bp{x) satisfies the Yang-Mills equations. The standard QCD La- 
grangian in this field looks like 



C = —Tr Gl + ]T (q f [ 7m (V m + igA,) - m f ] q f ) (2) 
G^{x) = V u Ap - V M A, + g[Ap(x),A v (x)\, Ap = iM£ 



where 



V M (x) = dp + iBp(x) = dp + iA 2 hbp U x u , 
[Vp(x), V u (x)} = -2iA 2 nbp V , n = n a t a . 

The part of the QCD Lagrangian which is responsible for meson hadroniza- 
tion can be written in the form 

C = (qS-'q) - ligD-'g) + g{qi 1 pt a q)g a p ) (3) 

where the quark field q(x) = qf aa ( x ) has indexes 

/ — flavor SUf(3), a — color SU C (3), a — spin 
and the gluon field g a n(x) has 

a — color SU C (3), n — vector. 
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The quark flavor spinor can be represented as 




SU f (3) : 



Quark propagator in the self-dual homogeneous gluon vacuum field is the 
solution of the equation (see details in [9]) 

(V(x) - m)S(x) = -S(x), V M (x) = <9 M + iA 2 hb fJiU x u . 



It has the form 

S±( P ) = 



du 

-e u ^ 



2A 2 



l-u 
1 + u 



(4) 



ip ± unj^bp) + 



m 



l-u 2 



The gluon propagator is obtained in [10]. We do not show it here because 
formula is quite cumbersome. 

We shall use the rough approximation of propagators preserving the main 
their features. The quark propagator having the Gaussian form and "zero 
mode" behavior for small quark mass ~ — is chosen in the form 

1 m. 



~ / \ ip 1± 1 



(5) 



UJ 



The gluon propagator is chosen in the form 

D%{x - x') = {ga»{x)g al »,{x')) = S aa 'S^D(x - x') 

where 



D{y) 



A 2 



(6) 



This rough choice of the quark and gluon propagators (5) and (6) is 
defined by the unique reason only: the Bethe-Salpeter equation can be solved 
analytically in this case and we get simple analytical formulas for the meson 
spectrum (see also [11]). 
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3 One-gluon exchange 

In one-gluon exchange the four-quark interaction Lagrangian is 

2 

W = ^ Jl dx 1 dx 2 (q(x 1 )i^t a q(x 1 ))D(x 1 -X2)(q(x 2 )i^t a q(x 2 )) (7) 

In order to go to colorless quark currents the Fierz transformations should 
be done 

• Color transformations 

4 1 

SU C (3) : (t a )j 1 j' 1 (t a )j 2 f 2 = g S jij 2 S 32j[ - ^ ta )hf^ ta )hf[ 

• Flavor transformations 

577/(3) : &hf'J>hf' 2 = £ S fi% s hfi + -^Whf'S^hfi 

• Dirac spin transformations 

Then for the four-quark interaction Lagrangian with pseudoscalar and 
vector colorless quark currents we get 

W P = dx l( ix 2 D( Xl - x 2 ) (8) 

• {^(9(^i)«759(^2))(fe)i75g(^i)) + ^(^i)^75A a g(x 2 ))(g(x 2 )i75A a ?(a; 1 ))| , 



W v = \-\ J/ dx l( lx 2 D( Xl - x 2 ) (9) 



Flavor structure of meson currents looks like 

= [cos 9{nn) + sin 8(ss)] 2 + [— sin 6{nn) + cos 6(ss)] 2 
+ (nr%)(nr%) + 2[(us)(su) + (ds)(sd)] 

where 

1 ( u \ / _ (uu) + (dd) 

n =W{d)> W = ^2 

The ideal mixing angle is 

did = arcsin —= = 35.3° 
v3 

Pseudoscalar nonet is defined as 

ty 3 = nr^n, K~ = us, (10) 

. uu + dd . n _ . uu 4~ dd n 
ri = cos v ■ = h sm t) ■ ss, n = — sin V ■ = h cos v ■ ss. 

Vector nonet is defined as 

fp = nT 3 n, K* = us, co = fin, (/) = ss. (11) 

4 Solution of the Bethe-Salpeter equation 

The solution of the Bethe-Salpeter equation is reduced to the following vari- 
ation problem (see [12]) 

Here 

nAp^ + yz) = Jdxe^C J Tr[r J S 1 ( X -^)r J S 2 (- X -^)~ 

<Ik e i ( fc ^+w))n J (fc,p) 



(2tt) 
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fl P (k,p) = Tr 



m2 48 



A 2 



fc 2 



M 2 

+ + 



A 2 1 



2A 2 8A 2 mim 2 



MM = -Tr 



il»Si k + 



V 



= -5 



_M 2 12 

^ A 2 



^71/52 ( /«• 
k 2 



M 2 

+ TTT7 + 



A 2 



2A 2 4A 2 m x m 2 



(1 - Wi^) 



A 4 



Tr : 4 • 3 • 2 = spin • color • duality 

To solve the variation problem (12) we choose the test function in 
form 

/ x aV , , r aV (2ir) 2 

U(y) = e~ a ^, (UU) = / dye- a ^~ = K ' 

J a' 

Then we get 

Let us define the vertex function 

i(kv) 4 a 

h«w — e o+i . 



= J dyV(y)e 



A(a + 1) 2 

The variation problem (12) is reduced to 

\j{M) = max(VnjV), J = P,V. 



9 



Let us consider 



with 



r dk 
J 



_ k A _ k z 2 
-e A 2 " A^'a+l 



2A 2 



A 4 (1 + a) 5 



(4tt) 2 (3 + a) 2 



1 + a 
3 + a 



We get 
F(A) = 



4 • 16 • 4£ 



max 



1 + a 



9-4 « \(l + a)(3 + a)J 

256 / a maa; (A) * 

3 \(l + a max {A)){3 + 



3 + a 
> 2 r 



1 + «m M (^) 
3 + a ma;E (A) 



+ A 



(16) 
F + Fi • A 



where 

F = 0.9188, Fi = 1.52 . 

The point a max (A) is defined by one of three roots of variation equation which 
is the algebraic equation of the third oder 



max 



(l + a)(3 + a) 



1 + a 
3 + a 



+ A 



-3 - 4a + a 3 + A(-9 - 3a + 3a 2 + a 3 ) = 
and can be found in the explicit form 



a max {A) ~ 1-7238 + 



0.3268 
A + 0.5646 



(17) 



Thus the eigenvalues of the Bethe-Salpeter equation can be written in the 
explicit analytical form. 

For pseudoscalar mesons one can get 



a 



-X P (M P ) 



* \ A 2 j m% (Ml 

e^F I ^+ + Hi 



7T 



8A 2 



(18) 
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where 



Hr 



H - Ai . l±A 

~~ 2 



A 



A 2 



*! m?, 

For vector mesons we get 



iZ _ = 4 ■ l±f COS 2 9 + 4 • ^ sin 2 # 

'/ mf. 2 mt 2 



-Ay(Afy) 



' e 7$-F \— + H v 



7T 



4 \4A 2 



where 



v 



H 



K* 



A 2 

m u m s 



[I - UJ U U S ) 



1 ^ = S-(i-^ 2 ) 



(19) 



5 Masses of pseudoscalar and vector mesons 

The experimental masses of pseudoscalar and vector mesons are (in Mev) 

M w = 140, M K = 495, M v = 547, M v > = 958, (20) 
M p = 770, M w = 782, M K * = 892, = 1019. 

We have four free parameters 

A, m u , m s , 9 (21) 

where A characterizes the confinement scale, m u = and m s are quark 
masses, 9 is the mixing angle. 
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The eigenvalues of the Bethe-Salpeter equation on the masses of our 
mesons look like 

-Ap(up) = 4 F ^— + H P ) , u P = ^— J (22) 



-\y(Vy) = ^lA e -i-lFl^-+H V ), V V = 



1 V 



2 



tt 4 V 4 V V A 

and these eigenvalues should satisfy the equations 

1 + X(v P ) = 0, 1 + X(v v ) = (23) 

We suppose that there exists the coupling constant a s (v), where a s (v) is 
a monotone decreasing function and in the points Vp and Vy this function 
satisfies 

a s (v P ) = irRp(v P ), a s (v v ) = ttR v (v v ) 

where 

Vp 

Rp M = tp ( vp \ V P = {*, K,V,r/}, (24) 
b {— + H P ) 

4 e~~ 

RvM = —j- — -y, V = {p, K*, u, 0}. 

b ^f + H v ) 

We consider that the experimental masses of pseudoscalar and vector 
mesons are fixed by (20). The problem is to find the parameters A, m u , m s , 9 
in such a way that a monotone decreasing function R(v) should be smooth 
as much as possible and satisfy (24). In our examples we have selected these 
parameters by eye and then have used the " Mathematica" program 

R(v) = Fit [{up, Rp(v P )}, {v v , R v (v v )}, {l,v,v 2 , ...,v n }, v] (25) 

for some n. 

Thus the coupling constant can be calculated 

a 9 (v) = nR(v) (26) 
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6 Effective coupling constant and meson cur- 
rents. 



The eigenvalues of the Bethe-Salpeter equation are the polarization operators 
of corresponding mesons. Let us consider the pseudoscalar mesons. We have 

U P (v) = X P (v) = R(v P )e%F Q + H^j (27) 

The kinetic term in the effective Lagrangian of pseudoscalar mesons looks 
in a vicinity of the mass shell 

(P[l - Xp(v)]P) - (P[l - Xp(v P ) - X' P (v P )(v - v P )]P) 

~ Zp (Mp-p>) 



{P[Zp(v P - v)]P) 



A 2 



(P[M P -p 2 }P), 



where the renormalization of the meson field is done 

A 



P(x) 



-.P(x) 



Here the constant of renormalization is 

d 



dv 



X P (v) 



X'p(v P ) 



V=Vp 



1 F'jvp) 
4 F(v P ) 



(28) 



The renormalization of the meson fields leads to the renormalization of 
the coupling constant and therefore to the renormalization of the vertexes 
which determine the binding of the meson fields with the quark currents. 

For the vertex the renormalization leads to 

gJ±=\v P (k) = -^=V P {k) = gpVp(k) 



gV P 



^Zp~A 



= 2ir\ 



la s (v P ) 



4a i 



irZp (ap + 1) 



-g A 1 ap + 1 



(29) 



where 



ap 



The interaction Lagrangian of the pseudoscalar meson P with quarks is 



2gp (- 



U{x) = 

with an appropriate matrix Ap. 



q{x)V P i-p) X P n 5 q(x) P(x) 



(30) 
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7 The decay constants fp 

The decay constants fp are defined by the matrix elements 

3?» = (0|J£(0)|P(Mp)) = V,- (31) 
The axial vector local currents with quantum numbers of pion and kaon are 



1 ,/w 1 



J ls( x ) = 7y| (u(x)^ 5 d(x)) , J^(x) = -j= (m(x)7 m 7 5 s(x)) 
The matrix element (31) is defined by the integral 

After simple calculations one can get 



2A /a,K) 12a* A 
/tt = — ~ 775- ^" e4 > ( 32 ) 

2A la s (v K ) 12a^ 1 r A A 

371 ' V kZ k ' n(2 + a K ) 2 ' 6 * ' 2 + ^J. 

8 Results and conclusion 

The numerical results of our calculations are shown on the Figures 1-3. They 
should be considered as preliminary qualitative estimations. The main con- 
clusions are 

• Spectrum of meson masses can be understood as result of (1) the ana- 
lytical confinement and (2) a monotone decreasing coupling constant. 

• Decay constants /„- and fx have reasonable values. 

• Different character of the behavior of a s (v) on Figures 1-3 means that 
mass spectrum of pseudoscalar and vector mesons can not fix parame- 
ters (21) uniquely. We should take into account other matrix elements 
(n — > 77, p — > 7T7T and so on). 
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Thus the simple form of quark and gluon propagators gives qualitative 
correct results, so that we can say that the structure of the quark-gluon 
interaction in the confinement region is guessed correctly. Besides we want 
to stress that relativistic scheme of hadonization was used and this scheme 
is very far from nonrelativistic potential picture. 



Another point is to connect the behavior of the coupling constant a s 



in the region M < 1000 Mev with known behavior in the region M > 
1000 Mev to compare with available results [6]. 

This work was supported by the grant RFBR 04-02-17370. 
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Figure 3: Pseudoscalar and vector mesons. Parameters A = 310, m l 
130, m s = 240, 9 = 64°. Decay constants f w = 131, f K = 172. 
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